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Abstract 



Ill-defined pinch singularities arising in a perturbative expansion in out of 
equilibrium quantum field theory have a natural analogue to standard scat- 
tering theory. We explicitly demonstrate that the occurrence of such terms 
is directly related to Fermi's golden rule known from elementary scattering 
theory and is thus of no mystery. We further argue that within the pro- 
cess of thermalization of a plasma one has to resum such contributions to 
all orders as the process itself is of non-perturbative nature. In this way the 
resummed propagators obtain a finite width. Within the Markov approxima- 
tion of kinetic theory the actual phase space distribution at a given time of 
the evolution enters explicitly. 
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Non-equilibrium many-body theory or quantum field theory has become a major topic of 
research for describing various transport processes in nuclear physics, in cosmological particle 
physics or more generally in quantum dissipative systems. A very powerful diagrammatic 
tool is given by the 'Schwinger-Keldysh' fl|-|3| or 'closed time path' (CTP) technique by 
means of non-equilibrium Green's functions for describing a quantum system also beyond 
thermal equilibrium Q. For an equilibrium situation this technique is equivalent to the real 
time description of finite temperature field theory 

Employing the diagrammatic CTP rules potential 'pinch singularities' might arise in 
strictly perturbative expressions. As an example we consider a scalar field theory. A typical 
contribution arising in a perturbative expansion takes e.g. the form 

Di ct (p, Po )i; (p,p )D™(p,p ) . (i) 



ret 




Here S describes some physical (perturbative) quantity (e.g. a self energy insertion); D 
and Dq v denote the free retarded and advanced propagator, respectively. As D^ et contains 
a pole at po = ±E P — ie and D$ v a pole at po = ±E P + ie the product of both in the above 
expression is ill-defined, if E (p, Po = E p = ^/m 1 + p 1 ) does not vanish onshell. Transforming 
such an expression back into a time representation, the contour has to pass between this 
pair of two infinitely close poles. 

It was observed and proven by Landsman and van Weert that such ill-defined terms 
cancel each other in each order in perturbation theory, if the system stays at thermal equi- 
librium [[J. Their arguments, however, rely solely on the KMS boundary conditions of the 
free propagators and self energy insertions, so that they do not apply for systems out of 
equilibrium. This severe problem arising for systems out of equilibrium was first raised by 
Altherr and Seibert ||. Indeed, it was speculated there that the CTP formalism might not 
be adequate for describing non-equilibrium systems at all. In a subsequent paper, Altherr 
H tried to 'cure' this problem by hand by introducing a finite width for the 'unperturbed' 
free CTP propagator D so that the expressions are at least well-defined in a mathematical 
sense. Within his modified perturbative approach, he also showed that seemingly higher 
order diagrams do contribute to a lower order in the coupling constant, as some of the 
higher order diagrams involving pinch terms will receive factors of the form 1/T n , n > 1 
reducing substantially the overall power in the coupling constant. In his particular case 
Altherr investigated the dynamically generated effective mass (the 'tadpole' contribution) 
within standard 4 — theory. (For the hard modes the onshell damping T is of the order of 
o(g A T).) Therefore he concluded that power counting arguments might in fact be much less 
trivial for systems out of equilibrium. We will come back to his observation below. 

In a recent work we have discussed in detail that modes or quasi-particles become 



thermally populated by a non-perturbative Langevin like interplay between noise and dis- 
sipative terms entering the non-equilibrium quantum transport equations. In the process 
of thermalization the full propagators necessarily must acquire some finite width (due to 
collisions or more generally due to damping). Plasmons behave as 'nonshell' modes p"T| . 



Strictly speaking, the evolution of a non-equilibrium system towards equilibrium is always 
non-perturbative. We will come back to this interpretation in more detail below. 

First, however, we will elaborate on the physical reason for the occurrence of pinch 
singularities in a strictly perturbative expansion, when an interacting system is prepared 
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with some non-equilibrium occupation of the particles. As a motivation we were inspired 
by the idea that in principle the Schwinger-Keldysh formalism is also adequate to describe 
simple scattering processes where e.g. only two initial particles are prepared at some fixed 
momentum states in the past. Hence, the perturbative scheme of the Schwinger-Keldysh 
formalism should give the same results as elementary scattering theory. 

To set the stage we start with some formulae and manipulations already presented in 
M. We follow the notation of [O. For simplicity we consider in the following a weakly 
interacting scalar 4 -theory. The initial state in the far past (assuming a homogeneous and 
stationary system) is prepared by specifying the momentum occupation number n(p) of the 
(initially non interacting) onshell particles. Note that this occupation number depends only 
on the three momentum p. (If specified with a thermal equilibrium distribution at some given 
temperature, n would be replaced by the onshell Bose distribution ub(E p = yjm? + p 2 ).) 
The occupation number n(p) enters the (free) propagator 

D<{p) = -2m sgn(po) S(p 2 - m 2 ) [Q(p )h(p) - 6(-p )(l + n(p))] (2) 
In addition, we note the form of the free retarded and advanced propagator: 

Dt /av (P) = ~ YT- • ( 3 ) 

To calculate perturbative corrections to the propagators we apply the Langreth-Wilkins rules 



12\ which are quite well-known within the context of the Schwinger-Keldysh formalism: 

D rct = D rct + /jrct^ret^ret + _ _ _ = . pret + ^ret ^ ^ 

D av = D™ + BgTgD w +...=: D* v + AD av , (5) 
D< = D< + DfZ^D^ + D T et Z<D™ + D<X™D™ +...=:£>< + A,D < , (6) 

where the dots denote multiple self energy insertions which we will not consider for the 
moment. Here the retarded and advanced self energies are given by the Fourier transforms 
of (cf. e.g. Q) 

S ret (x!,x 2 ) :=e(t!-* 2 ) [S>(x 1 ,x 2 )-S<(x 1 ,x 2 )] , (7) 
S av (x 1) x 2 ) :=Q(t 2 -ti) [S < (x 1 ,a; 2 )-S > (x 1 ,a; 2 )] . (8) 

The self energies S > and S < are related as S > (xi,x 2 ) = S < (x 2 ,Xi) in case of a scalar 
field theory. The self energy insertion So in a strictly perturbative expansion is given by 
a convolution of the initial free propagators. If the initial momentum distribution entering 
the propagator (|) is given by the Bose equilibrium distribution, the important relation 

E>(p) = e p °/ T S<(p), (9) 

holds, which is nothing but the KMS boundary condition. It is worth mentioning that our 
conventions are chosen such that iS < (p) is always real and non negative. In a transport 
theory (see below) it can be interpreted as the production rate for modes with the respective 
energy. As a characteristic example we discuss in the following the 'sunset' graph arising 
in scalar 4 -theory. This diagram is illustrated in fig. [TJ We choose this particular graph 
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as an example since the self energies Eg (p,Po — E p ) do not vanish onshell for thermal 
distributions (see e.g. pB|JT^ ]). This also holds for any non-equilibrium distribution n as long 
as the individual two-particle scattering contributions are kinematically allowed. Within 
finite temperature field theory the imaginary part of the self energy ('cut' diagram) taken 
onshell is connected to the scattering rate (as an illustration see fig. 0). On the other hand, 



there exist certain self energy insertions like the so-called hard thermal loop self energy [15 



or other one-loop diagrams |T^] which vanish on-shell due to simple kinematical constraints 
and thus do not cause any pinch problem. 

By inspecting (§-§[) more closely one finds that the perturbative corrections AZ} ret//av to 
the free retarded/advanced propagator are free of any pinch singularities as the emerging 
poles are all located at the same side of the contour. We note in passing that this also holds 
for multiple self energy insertions in (1,0) (see e.g. [Kj). In contrast, all three contributions 
to AD< are ill-defined. Using the identity 

vrsgn(po) S(p 2 - m 2 ) = ~ (D^( P ) - D™(p)) (10) 

together with (|]) we can further manipulate the three contributions of AD < by employing 
the Fourier transforms of the definitions ([7]) and (H). We find 

AD<(p) = AD< g (p) + AD< nch (p) (11) 

with a regular part, 

AD< g (p) = [e{p )n(p) - Q(- Po )(l+n(m 

x (Df (p)^(p)D*\p) - D»(p)I%'(p)D»(p)) , (12) 

and the part carrying the pinch singularities, 

A^plnchW = DT(P) [e(po) ((1 + n(p))^(p) - m*%(p)) 

+ Q(-p ) ((1 + n(p))^(p) - n(p)E<(p))] D^(p) . (13) 

The last expression is ill-defined, if the terms in the square brackets do not vanish onshell as 
already pointed out in M. The expression in the square brackets is familiar from standard 



kinetic theory (see e.g. |{3|,|T0"11 ) : Apart from a trivial factor one can interpret 

r*f(p) := [(1 + n(p))^o(p) ~ m^o(p)) (14) 

^J-jp po=E p 

as the net effective rate for the change of the occupation number per time. For an equilib- 
rium situation the occupation number is given by the Bose distribution and the self energy 
insertions fulfill the KMS condition @. Hence, for the equilibrium case the whole bracket 
exactly vanishes and no pinch singularities emerge. In contrast, this is not the case for a 
general non-equilibrium configuration ||. 

To shed first some light on the physical interpretation of this ill-defined expression one 
has to ask for observables which are affected by this singularity. Within standard scattering 
theory one would think about the probability for a particle of some initial momentum state 
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to be scattered into another momentum state. Therefore we ask, how the occupation number 
h has changed after a long time. The occupation number for the out-states can be readily 
extracted from D < by means of the formula (for a derivation see ||10|| ) 



n(p, t 



\ (° ut ) 



t (out) (out) 



Ep 

. 2 i " 



1 d d 



d d 



2E n dt dt' + 2^ dt dtf'iV 



d 3 x / d 3 y e ipx e- m {iD< (y, t; x, t')) 



t'=t 



i x i dpo (E p pi \ < 



(15) 
(16) 



When inserting (JTTJ) one finds by contour integration that AD< g only yields a finite con- 
tribution. The same holds true for the 8(— po)-term in (|13|) since the 'particle projector' 
+ Po) vanishes on the antiparticle mass shell. However, the O(po)-term of the 



'JL -|_ 



A. 

2E„ 



ill-defined expression AD< inch gives rise to the following infinite expression 
An(p) (out) = T eS (p) ■ 2n5(0) + finite contributions. 



(17) 



From analogy to the standard derivation of Fermi's golden rule in elementary quantum 
scattering theory one is immediately tempted to interpret this 6(0) singularity as the elapsed 
scattering time T — > oo. Indeed, this interpretation has very recently been conjectured by 
Niegawa in W7\, where he was also elaborating on the issue of pinch singularities in non- 
equilibrium quantum field theory. His major point, however, has been to interpret the 
infinite shift An(p) as a renormalization in the number density. We think, however, that 
this latter interpretation further obscures the problem instead of uncovering the physical 
processes which are at the bottom of the pinch problem. 

To demonstrate that the pinch singularities indeed appear as a result of Fermi's golden 
rule in scattering theory we now assume that the interaction is switched on at a time 
t = —T/2 and switched off at t = T/2, i.e. we replace 



S </>( Xl) X2 ) ^ f%i>{ Xli x % ) := 9(| - h) 8(f - h) ^ /> (x 1 ,x 2 ) G(h + |) G(t 2 + \ ) (18) 

and assume that the duration time T is large but finite. This procedure regulates the pinch 
singularity to a finite value. As a first step we again extract the pinch term from (|5p, now 
working in the representation of three-momentum and time: 

A^pinch^M') 

<#>0(1) <#> (2) rf P0(3) 



2vr 2tt 2tt 



x '(e(p 0(1 ))n(^-e(-p (i))(l + n(p)))SS v (^po(2)) + S<(p,p (2); 
- Ef (p,p (2)) (e(po( 3 ))n(p) - e(-po ( 3))(l + n{P)))[ 

T/2 T/2 
-T/2 ~T/2 



(19) 
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As 



T/2 

j dte itAp = A sin ^Ap) T ^°° 2vr5(Ap) 

-T/2 



(20) 



it becomes clear how the pinch singularity arises for T — > oo. Furthermore, if T is already 
sufficiently large, we are safely allowed to approximate p (i) ~ Po(2) ~ Po(3) within the square 
bracket in (|19|): 

(21) 



©(P0(2)) ((l+w(p))E^(p,p (2)) -^(p)S>(p,p (2), 

+ ©(-Po(2)) ((! + ™(r>) s o (P>Po(2)) -Hp)Zo(p,Po(2), 



We proceed by calculating An(p)( out ) by means of (0). For this we first take t, if > T/2, 
evaluate the po(i) - an d Po(3)- m t e g ra ti° n by standard complex contour integration and then 
insert the emerging expression into ([15|). It results in 



An(p)S ut) 



pinch 

dp (2) 
2ir 



1 d d i . d d . 
2E~ n dtdt' + 2^3t ~ dt'> 



©(Po(2)) ((1 + ™(p)) s o (P>Po(2)) -n(p)Ho(p } po(2) 

+ ©(-P0(2)) ft 1 + ™0?)) S O (P>f>0(2)) - «(p)E^(p,Po(2) 



X 



p L 



sin(|(p (2) + E p ))^ Ept sin(f (p (2) - 



X 



sin 



P0(2) + -Ep 

(P0(2) - #p)) 



P0(2) -p 



sm( 



T (p (2)+^p)) 



!''> e -iE p t 



P0{2) - E p 



P0(2) + E p 



I'D -iE p t' 



t'=t 



r e fr(p) 



dpo(2) 

2vr (p (2) 



^p) 2 



T 



sin 77(^0(2) = r cff (p) ■ t 



(22) 



which is valid for large but finite T. 

Thus we have demonstrated the bridge between the occurrence of pinch singularities 
within the context of the CTP formalism and Fermi's golden rule in elementary quantum 
scattering theory. The effective rate r e g is therefore analogous to the transition probability 
per unit time. Indeed one can easily understand in physical terms that one has to expect 
such a singularity in perturbation theory: Staying strictly within the first order contribu- 
tion the particles remain populated with the initially prepared non-equilibrium occupation 
number (since this quantity enters the free propagator (§)) and scatter for an infinitely long 
time. Therefore, the resulting shift An(p)( ont ^ should scale with r e g(p) • T with T e g(p) held 
fixed. We conclude that the occurrence of pinch singularities appearing in perturbative 
contributions within non-equilibrium quantum field theory is of no mystery, but actually it 
has to appear because of a very intuitive reason: the interaction time T becomes infinite. 
However, looking at a Boltzmann equation which describes the time evolution of the particle 
distribution function in the semiclassical regime (see (|30|) below) one realizes that the oc- 
cupation number does not stay constant during the dynamical evolution of the system, but 
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will be changed on a timescale of roughly 1/r. The quasi-particles are not really asymptotic 
states. 

Next, however, we will show how pinch singularities are formally cured by a resummation 
procedure. The onshell non-equilibrium effective rate T e g can be visualized as being the net 
result of collisions between the onshell particles. From standard thermal field theory one 
would thus expect that the propagators will become dressed and supplemented by a finite 
(collisional or more generally damping) width. This represents already a non-perturbative 
effect which only can be achieved by a resummation of Dyson- Schwinger type. As a first 
attempt (proposed by Baier et al. one might resum the full series of ([§-0) using the 

self energy S (recall that the latter is calculated from free propagators): 

D = D + D Q Z D Q + D Q Z D Q E Q D + . . . = D + D E D . (23) 

With the definitions T (p,p ) := ^[E^{p,p ) - Ejf&Po)] and Re£ := ReSff = Re£g v we 
end up with (cf. e.g. ||10|1 ) 



D ret = D r C t + jyfjtfjy* = (24) 

p 2 — m 2 — ReE + Wo^o 
D- = D? + fl-iyi,- = ^_ m ,_^_ w • (25) 

D< = Z^D" = (-* V _ m2 _^ o)2+?gr , . (26) 

Hence the resummation of the series (@) of ill-defined terms results in a well-defined expres- 
sion. The quantity 

appearing in fl26|) has to be interpreted as the 'occupation number' demanded by the self 
energy parts |1U|| . If the equilibrium KMS conditions (^) apply for the self energy part, 

then n,Y;(p,po) ^— -* n B (p ) becomes just the Bose distribution function. For a general 
non-equilibrium situation, however, this factor deviates from the Bose distribution. If the 
damping width is sufficiently small, i.e. f, S > , S < are proportional to some power in the 
(small) coupling constant g (e.g. ~ g A in case of the sunset graph depicted in fig. P the 
expression fl26|) results in 

D^YfiD™ ^4 -2m sgn(cj) 5(p 2 - m 2 ) lim n^(p,p ) . (28) 

When evaluating the occupation number for the out-states by means of (|I6|) one accordingly 
will get 

n(p)^ « n s (p,£ p ) (29) 

which is free of any pathological behavior. The astonishing thing to note at this point is 
that in fact the (initial) non-equilibrium distribution n has been substituted by and, 
therefore, does not show up explicitly. So the question is, how ft enters? 
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Calculating Eo on a purely perturbative level the initial occupation number h enters via 
the free propagator ([|). This however cannot be the whole truth in a dynamically evolving 
system. It is important to make sure that such a system is prepared at some finite initial time 
t . (If t would be taken as t — > — oo the system would already have reached equilibrium 
long time ago.) Bedaque [|19[] already has noted that pinch singularities are in fact an artifact 
of the boundaries chosen at to — > — oo. Time reversal symmetry is explicitly broken, so that 
the propagators in principle have to depend on both time arguments explicitly before the 
system has reached a final equilibrium configuration. Therefore the use of Fourier techniques 
(which in fact has led to the pinch singularities in flT3p) is highly dubious. The initial out 
of equilibrium distribution n(t ) cannot stay constant during the evolution of the system as 
it has to evolve towards the Bose distribution. Hence there must exist contributions which 
attribute to the temporal change of the distribution function. As long as the system is 
not in equilibrium (on a time scale of roughly l/T Q (p, E p )), the propagator thus cannot be 
stationary. In addition, the self energy parts S < and S > do also evolve with time. Hence 
they should depend on the evolving distribution function and not persistently on the initial 
one, n, which enters So in (|23|). Thus the resummation of fl23|) does not cover all relevant 
contributions. Speaking more technically, the self energy operators must also be evaluated 
consistently by the fully dressed and temporally evolving one-particle propagators. 

The solution to these demands is, of course, the description of the system by means of 
appropriate (quantum) transport equations |20|,|4],|10H . Graphically this is illustrated in fig. || 
In addition to the sunset diagram we have also included the mean field or Hartree diagram 
there which in a perturbative scheme is the one which arises first. (It would, however, not 
result in a pinch singularity so that we had discarded it in our previous discussion.) The 
difference to the resummation of ( p3|) is the fact that the propagators entering into the 
self-energy operators are now also the fully dressed ones. Such a skeleton expansion of the 
self energies with including the dressed propagators in the resummation is also familiar in 
standard quantum many-body theory for strongly interacting systems pifl . 

Unfortunately, the full quantum transport equations are generally hard to solve and 
thus are not so much of practical use. Yet one need not be that pessimistic. If the coupling 
is weak, i.e. the damping width is sufficiently small compared to the quasi-particle energy 
(which one typically assumes for many applications) one can take the Markov approximation 
to obtain standard kinetic equations (for a derivation see e.g. [|20|,fiyTO| ) . For the situation 



illustrated in fig. |3] one gets the standard form fTO 



(E p dt - pdg - d 3 m(x, t)d^) f(x, t; p) (30) 
= 2 [^(X't'iPiEp) (f&P,t) + !) -«£>(£, E p )f(x,p,t)] 

Here / denotes the semi-classical non-equilibrium phase-space distribution of quasi-particles. 
m(x, t) denotes the sum of the bare and the dynamical (space time dependent) mass gen- 
erated by the Hartree term. Within the spirit of kinetic theory one easily realizes that the 
result obtained in (p2|) simply states that the change in the occupation number per time T 
is nothing but the collision rate. Within this Markovian approximation the fully dressed 



propagators are given by |T0] 



D TCt (x,t;p) « - (31) 

pz _ m ^( Xl t) — ReL(x, t; p) + tp l [x, t; p) 
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v ' " ' p 2 — m 2 (x, t) — Re£(x, t; p) — ip T(x, t; p) ' 
D<(x,t;p) « ("20 (p2 _ m2( . t) _ ^° eS ( ( ^ ;.^ ))2 - ^ r2(f ; t . p) /(*. *! P) • ( 33 ) 



In particular we emphasize that in (|33|) the instantaneous non-equilibrium phase space dis- 
tribution function f(t) enters and not the initial one, n. The dynamically generated mass 
as well as the collisional self energy contribution £ can thus be evaluated with these propa- 
gators. (In kinetic theory one usually takes the propagators in their quasi- free limit (ReS, 
r — > 0), albeit instantaneous.) Higher order terms leading to the pinch singularities are 
explicitly resummed and lead now to finite and very transparent results. 

One can now easily understand the observations made by Altherr ||. He has found, 
starting from some non-equilibrium distribution n, that higher order diagrams contribute to 
the same order in the coupling constant as the lowest order one. Indeed, in his investigation, 
the particular higher order diagrams where nothing but the perturbative contributions of the 
series in (Q) for the dressed or resummed one-particle propagator D K . The only difference 
is that he has employed a 'free' propagator modified by some finite width in order that 
each of the terms in the series (|J) becomes well defined. The reason for the higher order 
diagrams to contribute to the same order is that the initial out-of-equilibrium distribution n 
cannot stay constant during the evolution of the system as it has to evolve towards the Bose 
distribution. If h — n% is of order o(l), it is obvious that there must exist contributions which 
perturbatively attribute to the temporal change of the distribution function and contribute 
to the same order o(l). In fact, in our prescription ([33|), n has simply be substituted by the 
actual phase space distribution /. Then calculating e.g. the tadpole diagram, as discussed 
in the particular case of [[J, one has to stay within lowest order in the skeleton expansion, 
but with the fully dressed propagator. 

In summary, we have shown in simple physical terms why so called pinch singularities 
do (and have to) appear in the perturbative evaluation of higher order diagrams within the 
CTP description of non-equilibrium quantum field theory. They are simply connected to the 
standard divergence in elementary scattering theory. The occurrence of pinch singularities 
signals the occurrence of (onshell) damping or dissipation. This necessitates in the descrip- 
tion of the evolution of the system by means of non-perturbative transport equations. In the 
weak coupling regime this corresponds to standard kinetic theory. In this case we have given 
a prescription of how the dressed propagators can be approximated in a very transparent 



form. Technically, pinch singularities appear due to a misuse of Fourier techniques p9 
From a physical point of view, scattering processes which change the occupation number 
give rise to pinch singularities, if these processes go on for infinitely long time. However, 
exactly these processes drive the system towards thermal equilibrium within a finite time 
characterized by the inverse damping rate. In equilibrium the occupation number stays 
constant and no pinch singularities can appear. 
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FIGURES 



FIG. 1. Lowest order self energy term in 4 -theory which contributes to the pinch problem 
(sunset diagram). 




FIG. 2. Imaginary part of the sunset diagram which can be identified with a scattering ampli- 
tude. 




FIG. 3. Dyson-Schwinger equation with fully dressed propagators (skeleton expansion). 
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